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Abstract 

QCD lattice simulations determine hadron masses as functions of the 
quark masses. From the gradients of these masses and using the Feynman- 
Hellmann theorem the hadron sigma terms can then be determined. We 
use here a novel approach of keeping the singlet quark mass constant in our 
simulations which upon using an SU (3) flavour symmetry breaking expan- 
sion gives highly constrained (i.e. few parameter) fits for hadron masses 
in a multiplet. This is a highly advantageous procedure for determining 
the hadron mass gradient as it avoids the use of delicate chiral perturba- 
tion theory. We illustrate the procedure here by estimating the light and 
strange sigma terms for the baryon octet. 
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1 Introduction 



Hadron sigma terms, ai^^ are definecQ as that part of the mass of the hadron 
(for example the nucleon) coming from the vacuum connected expectation value 
of the up (m) down (d) and strange (s) quark mass terms in the QCD Hamiltonian, 

= mf{H\{uu + ddY\H) , (Tp) = mf(/7|(ss)^|if) , (1) 



where we have taken the u and d quarks to be mass degenerate, mu = rud = 
mi. (The superscript ^ denotes a renormalised quantity.) Other contributions 
to the hadron mass come from the chromo-electric and chromo-magnetic gluon 
pieces and the kinetic energies of the quarks, [2]. Sigma terms are interesting 
because they are sensitive to chiral symmetry breaking effects. Experimentally 
the value for a^^"^ has been deduced from low energy vr-A^ scattering. A delicate 
extrapolation to the chiral limit pj gives a result for the isospin even amplitude of 
(^■kN I fl (with a-j^N = ^l^^)! from which the sigma term may be found. The precise 
value obtained this way has been under discussion for many years. However 
within the limits of our lattice calculation, this will not concern us here and for 
orientation we shall just quote a range of results from earlier analyses of [3l H] 
of 45(8) MeV while a later dispersion analysis [5j suggested a much higher value 
64(7) MeV. An estimation using heavy baryon chiral perturbation theory gave 
45 MeV, [6]. A more recent estimate gave 59(17) MeV, [7j. Even less is known 
about the nucleon strange sigma term. Eq. ([1]) is usually written (in particular 
for the nucleon) as 

^(N) ^ mf{N\{uu + dd-2ssr\N) ^ 2{N\issr\N) 

' ' ^ {N\{uu + dd)^\N) ' ^ ' 

(i.e. we consider a^^^ and ?/(^)^ rather than a^^^ and ai^'^). The simplest cal- 
culation, e.g. [1] (which we will discuss in more detail later) uses first order in 
SU (3) flavour symmetry (octet) breaking to give 

(;v) ^ mf Mh + Ms - 2Mm 26 

mf-mf 1 - ^ 1 - ' 

and 

) = ~ 325 y MeV , (4) 

where mf/mf is the ratio of the strange to light quark masses, which using the 
leading order PCAC formula for this ratio gives 

/mf = (2M|, - M'J/M'^ ~ 25 . (5) 



^Or more accurately as the matrix element of the double commutator of the Hamiltonian 
with two axial charges. However this is equivalent to the definition given in eq. (P), see for 
example [1]. 
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The Zweig rule, (A^|(ss)-"|A^) ~ would then give 

al^^ ~ 26 MeV , af^ ~ MeV , (6) 

while any non-zero strangeness content, > would increase this value of 

crl^\ ai^^ (and indeed, due to the large coefficient, ai^'^ quite rapidly). 

Determination of the strange sigma term (and in particular ?/(^)^) is impor- 
tant in constraining the cross section for the detection of dark matter. WIMPs 
would be scattered off nuclei by the exchange of scalar particles, such as the Stan- 
dard Model Higgs particle, which will interact more strongly with heavier quark 
flavours. This coupling can be parameterised in terms of the fractional contribu- 
tion of a quark flavour q to the nucleon's mass Mjy, /j-^ = {N\{qq)"'\N) / . 
While the contributions of the charm and heavier flavours approach a constant 
that is proportional to the gluonic contribution /t^, there is a strong dependence 
of the cross section on the value of fx^ , see e.g. [S], [H] and references therein. 

Computing the sigma terms from lattice QCD has a long history from initial 
quenched simulations to 2 flavour and more recently 2 + 1 flavour simulations, 
e.g. [T0l[IIl[T2l[l3l[Tl[T5l[T6l[I7l[T8l[l9ll^ with a status report being given 
in [22]. In general more recent results tend to give a lower ai'^^ term than earlier 
determinations. 

In this article, we shall investigate this simple picture as described in eqs. ([3]), 
and in particular test the linearity assumption of SU{3) flavour symmetry 
breaking. 



2 Flavour symmetry expansions 

Lattice simulations start at some point in the {mf, mf) plane and then approach 
the physical point {mf*,mf*) along some path. (In future we shall denote the 
physical point with a *.) As we shall be considering flavour symmetry breaking 
then we shall start here at a point on the flavour symmetric line mf = mf and 
then consider the path keeping the average quark mass constant, m = const.. 
The SU{3) flavour group (and quark permutation symmetry) then restricts the 
quark mass polynomials that are allowed, [23], giving for the baryon octet 

Mh = Mo{m) + CHSmi + 0{6mf), (7) 

with 



ch = 



3Ai 


H 


= N 




H 


= A 


-3A2 


H 


= S 



(8) 



[ -3(A - A2) H = E 

where 

5mi = mi — rn , m = \{2mi + m^) , (9) 
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and Ai and A2 are unknown coefficients. So to linear order in the quark mass, 
we only have two unknowns (rather than four). A similar situation also holds for 
the pseudoscalar and vector octets (one unknown) and baryon decuplet (also one 
unknown). These functions highly constrain the numerical fits. (At 0{5mf) only 
the baryon decuplet has a further constraint.) 

Permutation invariant functions of the masses Xs-, (or 'centre of mass' of the 
multiplet) can be defined which have no linear dependence on the quark mass. 
For example for the baryon octet we have 

Xn = 1{Mn + Ms + Mh) = Mo(m) + 0{6m^) . (10) 

(The corresponding result for the pseudoscalar octet is given later in eq. fl29|) .) 

Furthermore expanding about a specific fixed point, mi = rris = friQ on the 
flavour symmetric line and allowing m to vary, we then have 

Mo(m) = Mo(mo) + Mo(mo)(m-mo) + O((m-mo)^). (11) 

We will see that Ai, A2 give all the non-singlet hyperon sigma terms and M'(mo) 
the singlet terms. 

As an example of the quark mass expansion from a point on the flavour 
symmetric line in Fig. [T]we plot the baryon octet Mh/Xn for H = N, A, H, S 
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Figure 1: Mh/Xn {H = iV, A, E, S) against Ml/Xl for an initial point ("sym. pt.") 
on the flavour symmetric line given by = 0.12090, left panel, and kq = 0.12092, 
right panel. The 32^ x 64 lattices are filled circles, while the 24^ x 48 lattices are open 
triangles. Also shown is the combined fit of eq. (j33p (the dashed lines) to the 32^ x 64 
lattice data. The fit results are the open circles, while the experimental points are the 
(red) stars. I and s denote the light and strange quark content of the hadron. 

against Ml/X"^ together with a linear fit, eq. (JTj) and implicitly eq. ( I29l) using 
2 + 1 0(a) improved clover fermions at /3 = 5.50, [24] using two starting values 
for the quark mass on the flavour symmetric line, namely kq = 0.12090, 0.12092. 
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All the points have been arranged in the simulation to have constant m. We 
see that a linear fit provides a good description of the numerical data from the 
symmetric point (where M^r ~ X* = 410.9 MeV) down to the physical pion mass. 

In a little more detail, the bare quark masses are defined as 

arug = - ( j , with q = l,s,0, (12) 

(with the index g = denoting the common quark along the flavour symmetric 
line) and where vanishing of the quark mass along the SU{3) flavour symmetric 
line determines kq-c- Keeping m = constant = mo gives 

= _3 _2_ • (13) 

Ko ft; 

So once we decide on a ki this then determines Kg. Note that kq-c drops out of 
eq. f lT3|) . so we do not need its explicit value. These initial Kq values chosen here, 
namely kq = 0.12090 and 0.12092 are close to the path that leads to the physical 
point (kq = 0.12092 being slightly closer). (This is discussed in more detail in 
[23], which also contains numerical tables and phenomenological values for the 
hadron masses. Results not included there are given in Appendix O) This path 
is also illustrated later in section 14. 3[ Fig. |H Although finite size effects tend to 
cancel in ratios of quantities from the same multiplet, we nevertheless fit just to 
the results from the 32'^ x 64 lattices (filled circles) using the linear fit of eq. (I7j). 
Finally note that we also have a similar flavour expansion for the pseudoscalar 
octet as for the baryon octet, as will be discussed in section H73l 



3 (Hyperon) scalar matrix elements 

Scalar matrix elements can be determined from the gradient of the hadron mass 
(with respect to the quark mass) by using the Feynman-Hellman theorem which 
is true for both bare and renormalised quantities. So if we take the derivative 
with respect to the bare quark mass we get the bare qq matrix element, 

^ = muu + dd)\H), ^ = {H\ss\H), (14) 

while if we take the derivative with respect to the renormalised quark mass we get 
the renormalised matrix element. In the left panel of Fig. [21 we show the nucleon 
masses (green diamonds) and the flavour symmetric nucleon masses (maroon 
squares) against l/ni, I/kq respectively (from eq. (fT2l) these are proportional to 
the bare quark mass). From the Feynman-Hellmann theorem, the slope of the 
masses (maroon squares) gives the total ^q=u ^ s(A^|gg|A^), while the slope of the 
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Figure 2: The left panel shows the nucleon mass, aMjsf, versus l/m, (for the 
fn = const, points, green diamonds with hq = 0.12090) and versus 1/kq (for the 
flavour symmetric points, "sym. pts.", maroon squares). The common flavour sym- 
metric points are denoted by red circles. The 24'^ x 48 volume results are open symbols 
together with a dashed line for the (linear) fit, while the 32^ x 64 volume results are 
filled symbols and solid lines. Similarly the right panel shows the nucleon mass aMjy, 
versus (aM^)^ (same notation as for the left panel). 



masses (green diamonds) gives the valence contributioiil. The difference between 
the two contributions gives the disconnected contribution. Because here all three 
quark masses are equal, the disconnected contribution for all three quarks will 
be the same. The two slopes thus give the estimates 

E,{NmN)con 4.0 
E,(iV|gg|iV) 9.7 

for bare lattice quantities. 

To look at renormalised matrix elements, we need a plot against the renor- 
malised mass, (aM^r)^ (as in leading order PCAC, is proportional to the 
renormalised quark mass, eq. (13T]) ). This is shown in the right panel of Fig. [2l 
The slopes are now much closer to each other. We now find the estimates 

V (iV|(gg)"|iV)con 3.2 



E,{Nmr\N) 4.3 

(A^|(ss)«|Ar) 1 /4.3 -3.2 



Ea{N\{qqr\N) 3 V 4.3 



0.085, (16) 



^Eq. d?]) can be extended to the 'partially quenched' case, [23], where the sea quark masses 
remain constrained by to = const, but tlie valence quark masses /i;, are unconstrained. 
Defining 6fiq = iiq — rn then for the nucleon, the leading change is particularly simple, CNdrrii — >■ 
cnS^i- For the other members of the octet, A, E, S, both Spu, 5^s occur, |23| . 
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for renormalised lattice quantities, giving ?/(^)^ ~ 2 x 0.085/(1 — 0.085) ~ 0.19. 
So although for bare matrix elements, there is a significant strange quark content 
this is reduced in the renormalised matrix element. 

We shall now try to make these considerations a little more quantitive. 



4 (Hyperon) a equations 
4.1 Renormalisation 

For Wilson (clover) fermions under renormalisation the singlet and non-singlet 
pieces of the quark mass renormalise differently [251 126] . We have 

= Z"' [rrig + az\{1mi + m,)] , az = — . (17) 

In the action the term '^qTri^q = ^ renormahsation group in- 

variant or RGI quantity. Upon writing this in a matrix form and inverting gives 



qq — \{uu + dd + ss) 



so for az ^ ^ then there is always mixing between bare operators. 

As an example of where this manifests itself, the relation between the bare, 
y'^^\ and renormalised y^^^"-^ cf. eq. (|2]), is then given by 

^ l + \az{l-y(^)) ' ^ ^ 

so we see that y^-^)^ ^ y^^^ for clover fermions. Additionally, since > and 
y^^^ > we find that < y^^\ i.e. is reduced. 

Useful quark combinations are the octet and singlet combinations, namely 

{uu + IdY - 2{'ssY = [{uu + dd) - 2(ss)] , 

Zj 

{uu + ~ddY + i^ssY = zNsn+ az) ^^^"^ ^ ^"^^ ^ ^^^^^ ' 



Furthermore, using the Feynman-Hellman theorem, eq. flMj) and with the hadron 
flavour expansion, eq. (171) together with eq. (tH]) gives 

{H\{uu + ddr-2{ssr\H) = (21) 

{H\{uu + ddr + {ssr\H) = ^T^. (22) 

Zj I -\- OlZ 

Eq. (12 ip . the equation for the matrix element of an octet operator, only involves 
cu (the hadron mass expansion keeping the singlet quark mass constant), while 
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eq. ( 122|) . the matrix element of a singlet operator, only involves Mq (occuring when 
changing the singlet quark mass). Eq. fl2Tl) also leads to eq. ([3]) as discussed in 
the introductioiij. 

Finally note that the quantities 

{tUs -mi){H\{uu + dd) -2ss\H) , {2mi + ms){H\{uu + dd) + ss\H) , (23) 

are RGI, all Z factors cancel when they are renormalised. Linear combinations of 
these two quantities are also RGI in particular the combination used previously 
of al^^ + as = ^^mq{H\qq\H). However, a[^^ and ai^^ considered separately 
are not RGI, see eqs. flT7|) . f fTSj) . The renormalised quantities are mixtures of 



the two lattice quantities, and is needed to relate lattice values to continuum 
values. Refering back to Fig. |2] we see that the bare lattice strange sigma term 
is much larger that the renormalised strange sigma term, due to a cancellation 
between the two terms in eq. flTSl) . 

4.2 a equations 

Multiplying the renormalised quark mass, eq. ( 1T7|) . together with eqs. (12^ 
(or more generally with eq. ( ITSll ) we can find RGI combinations (i.e. a form where 
the renormalisation constant Z'^^ cancels). In particular we find 

a;^)-2mf) = j^il + az)moCH (24) 

af) + raf) = ^^moM^(mo) , (25) 

where r is the ratio of quark masses 



r^^. (26) 



Thus we have to find the (fixed) coefficients (1 + az)'moCH, moMQ^mo) . We then 
determine the physical values of the sigma terms by extrapolating to the point 
where the quark mass ratio takes its physical value, i.e. r = r*. 

We observe that we have two simultaneous equations, which can be easily 
solved to giv^ 

o-J^^ = [(1 + az)"^oC/f + 2moMo(mo)] 

) = [-(1 + (^z)moCH + moM^(mo)] . (27) 



^The RHS of eq. ^ can be re-written as cn/Z'^^ = SAi/Z^-^. Together with Ms + My. - 
2Mn = —'dAiSmi = 3^1 (mf — mf-)/Z^^ this gives eq. ([3]). An ahernative mass combination 
that also picks out the Ai coefficient is Mh — Ma = —iAiSrai. 

^This leads to relations between the various sigma terms, which we list in Appendix \K\ and 
where we also argue that they are always approximately true. 
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We see that the smallness of cr| in comparison to as is certainly guaranteed by 

( H) 

the presence of an additional r in its numerator. As o"s > we must also have 
MQ(mo) > (1 + az) "oiaxcH- These coefficients are also sufficient to determine 
y{H)R^ as can be seen either directly from eq. f l27|) or from eq. f l22|) . 

y{H)R ^ 2 -(1 + (^z)moCH + moM^(mo) ^^^^ 
(1 + az)moCH + 2moMo(mo) 

Again, as seen in section |31 y^^^'^ only depends on gradients and not on the 
physical point. 

It is now convenient to normalise the coefficients by Xn so we now need to 
find the coefficients (1 + az)TnoCH / XN^mo) and moMQ(mo)/XAr(mo). 



4.3 Determination of the coefficients 

The hint for determining the coefficients from our lattice data is given in section |3l 
where we consider gradients with respect to a renormalised or physical quantity 
- here taken as the pion mass. As in eq. ([7]) we also have a similar expansion for 
the pseudoscalar octet, 



Ml = Ml^ + 2a5mi + 0{5ml 



(29) 



(together with M|- = M^^ - a5mi + 0{5mf), = M^^ - Aa5mi + 0(5mf)). 
This gives a good representation of the data as can be seen from Fig. 12 of [23] . 
Analogously to eq. ffTOj) we can define a flavour singlet quantity 



Xl = \(2Ml + Ml) = M^ + Oi6m^) . 
However, as well as eq. ((71), we have the additional constraint from PC AC 



(30) 



Ml 



2B^mf , 



(31) 



(together with M^ = B^{mf + mf ), M^^ 



2i?Q ) which implies that 



Ott 



2q;(1 + az)m . 



a 



(32) 



If we now consider an expansion in the (physical) pion mass then eliminating 5mi 
between eq. ([7]) and eq. ( 129|) gives 



M 



H 



X 



N 



'I + az)mo 



Ch 
X 



N 



+ 



1 + az)mo 



Ch 
X 



N 



Ml 



(33) 



from the point on the symmetric line mo = m. Thus if we plot Mh/Xj^ versus 
Ml/Xl (holding the singlet quark mass, m constant) then the gradient imme- 
diately yields (1 + az)moCH/XN. The only assumption is that the 'fan' plot 
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splittings remain linear in Smi down to the physical point. In Fig. [T]we show this 
plot giving the results 

{l + az)mo—^= 0.1899(55), 0.2066(68), 
Xn 
3A 

(1 + az)mo—^ = 0.03942(314) , 0.04164(431) , (34) 
Xn 

for Ko = 0.12090, 0.12092 respectively. 

Alternatively on the flavour symmetric line, mi = m (i.e. 6mi = 0), so varying 
m from a point mo gives 

M2(m) = Mo\(m) = M^^{mo) + M^l{mo){m - rrio) 

= 2a{l + az)[mo + {m - mo)] , (35) 

which gives M^^lmQ) = 2a{l + az)- So now eliminating (m — mg) between 
eqs. ffTT]) . (|35D gives 

XNjm) / 

XN{mo) V 

Again in a plot of Xj^ijn) /Xj^^mQ) versus X^(m)/X^(mo) the gradient immedi- 
ately gives the required ratio moMQ(mo)/Xjv(mo). We have also replaced Mjy by 
X]sr and by (which allows us to use all the 32'^ x 64 data available for a 
particular m). In Fig. |3]we plot XN(fn) / Xj\f{mo) versus X^{m) / X^{mo) . From 

1.10 
1.05 
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£ 

^ 0.95 
lE 

z 
X 

0.90 
0.85 
0.80 
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X/( m)/X/(m„) 



moM(^(mo) 

XAr(mo) 



+ 



moMo(mo) 

XAr(mo) 



xm 

Xl{m,) 



(36) 




Figure 3: Xjv(m)/Xjv(mo) versus X^(m)/X^(mo) along the flavour symmetric line, 
together with the linear fit from eq. (j36p . 



eq. fl36l) this gives 

moM^(mo) 
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Finally the quark mass ratio, r, must be estimated. In Fig. |l]we plot (2M^ — 




0.0 I ' ^ 
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Figure 4: (2M|, - M^)/X% versus M^/X% for kq = 0.12090 (left panel) and kq = 
0.12092 (right panel). The 32^ x 64 volume resuhs are given by the filled symbols, 
while the 24^ x 48 volume results are shown using empty triangles. The fit is given in 
eq. (f38l) . Experimental points are denoted by (red) stars. 

Ml)/Xj^ versus M^/Xj^. From eq. ([29]) we have 

2Ml-^_ Ml Ml 

As in section |2l we see that for constant m the data points lie on a straight line 
(i.e. there is an absence of significant non-linearity). Furthermore the gradient is 
fixed at —2. (Indeed leaving the gradient as a fit parameter for the = 0.12090 
confirms that this gradient is very close to —2.) Together with PCAC, eq. fl3T|) 
this gives the x-axis is proportional to mf while the |/-axis is proportional to 
and thus the ratio gives r. Taking our physical scale to be defined from M^/X^|* 
(i.e. from the x-axes of Fig. H]) gives 

27.28(16) Ko = 0.12090 , . 

26.23(24) Ko = 0.12092 ' ^ ' 



mf 



4.4 Curvature effects 

What can we say about corrections to the linear terms? The simple linear fit 
describes the data well, from the symmetric point to our lightest pion mass, both 
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along the m = const, line and the flavour symmetric line. To see qualitatively 
the possible influence of curvature we now compare linear flts with quadratic flts. 
These will be used to estimate possible systematic effects. We briefly discuss 
these effects here. 

In Fig. Owe compare the results of a quadratic flt and a linear flt, both for the 
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Figure 5: Left panel: Mh/^n for H = N , A, 12, r, against aSnii for initial point on 
the flavour symmetric line given by kq = 0.12090 together with the previous linear fit 
(dashed lines) and quadratic fit (solid lines). Other notation as in Fig. [H Right panel: 
Xj^{rn)/XN{mo) versus X^(m) / X^{mo) along the flavour symmetric line, together 
with a linear fit from eq. (f36]l (dashed line) and a quadratic fit (solid line). 



baryon mass fan plot and for Xiy(fn)/Xj^{mo). In the left panel of the figure, we 
consider the baryon mass fan plot. The quadratic fit here uses all the data, [23] . 
on both lattice sizes (in cases where results for two lattice sizes are available, 
we used the larger lattice size only). The curvature terms here are small and 
statistically compatible with zero. 

The right panel of the figure shows a quadratic fit to the results along the 
symmetric line. The curvature here is dominated by the large error of the lightest 
point (which has a low statistic). Thus we shall regard this fit as only giving an 
estimation of the possible systematic error. 

The results in the next section include systematic error estimates from both 
these curvature sources combined in quadrature. In Appendix [B] we give some 
more details. 



5 Results 

We can now numerically determine y^^^^ and cr;^"*, (ri^\ 

We start with ?/(^)^. From eq. together with eqs. (EID, (EZD and eq. dS]) 

gives the results in Table [TJ The first error is the linear fit error (in this case 
dominated by the error in eq. (!37|) ). while the second error indicates possible 
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N A S S 

/to = 0.12090 

'^(^* 0.22(9)(15) 0.80(14)(28) 1.23(20)(41) 2.14(38)(64) 
a{^^*[MeV] 29(3)(4) 23(3)(4) 20(3)(4) 16(3)(5) 
af'^*[MeV] 89(34)(59) 250(34)(68) 334(34)(68) 453(34)(58) 
kq = 0.12092 

Ip^* 0.18(9)(15) 0.79(14)(28) 1.25(20)(42) 2.30(42)(68) 

cT;^)*[MeV] 31(3)(4) 24(3)(4) 21(3)(4) 16(3)(4) 
cTp^*[MeV] 71 (34) (59) 247(34) (69) 336(34) (69) 468(35) (59) 



Table 1: Results for the baryon octet for y^^^^*, aj^^*, ar'* with H = N, A,T.,E 
for Ko = 0.12090, 0.12092. 

effects from liigfier order terms, as discussed in section 14.41 We see tliat tliere 
is an order of magnitude increase in tfie fraction of {H\(ss)'^\H) compared to 
{H\{uu + dd)^\H) as we increase the strangeness content of the baryon from the 
nucleon (no valence strange quarks) to the S (two valence strange quarks). 

Turning to the sigma terms themselves, from eq. fl2^ we can find an indication 
of the magnitude of a\^^ as approximately (with Xjy = 1.1501 GeV), 

{N)* 

^{N) * _ [22 ~ 25] + MeV > [22 ~ 25] MeV , (40) 

(for Ko = 0.12090, 0.12092 respectively). The last inequality follows as obviously 
^(Af)* ^ Indeed this shows that a non-zero ai'^^* > can only add a few MeV 
to this result. 

The results for o'l^''* and ai^^* are also given in Table [H (Again the first 
error is the statistical error, while the second systematic error is due to possible 
quadratic effects.) While the data for kq = 0.12090 is more complete than for 
Ko = 0.12092 (cf. the plots in Fig. [1]) and demonstrates linear behaviour, as the 
path starting at ko = 0.12092 is closer to the physical point (cf. Fig. Hj) we shall 
use these values as our final values. These results are illustrated in Fig. [6] for 
y{H)R* where H = N, A, S, S. 

By varying r in eq. fl27j] |^ , we plot in Fig. [7| aj^^ and ai^^ for the baryon 
octet, H = N, A, and S from the symmetric point (vertical dashed line at 
X = 1) to the physical point (left vertical dashed line), a^^^ is rapidly decreasing 

^Using, for example, the resuhs from the left panel of Fig. |4l r may be re-written as 

3 - 2(M2/X2) • 
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Figure 6: y(^)«* for = iV, A, S, H using the results from Table □ for kq = 0.12092. 
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Figure 7: (decreasing from the symmetric point x = 1) and al (increasing) for 
H = N , A, T,, for kq = 0.12092. The physical and symmetric lines are denoted by 
vertical dashed lines. 
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while (t1^^ is increasing as we decrease the quark mass. Also, as expected a}^'' is 
largest for the nucleon, A^, while ai'^^ is the smallest. Finally in Fig. [8] we plot 
(tI^^* , cri^^* against H = N, A, J] and S, again using Table [H 




Figure 8: a} '* and ar'* for H = N, A,T,,E at the physical point for kq = 0.12092. 



6 Conclusions 

Keeping the average quark mass constant gives very linear 'fan' plots from the 
flavour symmetric point down to the physical point. This implies that an expan- 
sion in the quark mass from the flavour symmetric point will give information 
about the physical point. In this article we have applied this to estimating the 
sigma terms (both light and strange) of the nucleon octet. There has been no 
use of a chiral perturbation expansion (indeed this is an opposite expansion to 
the one used here, expanding about zero quark mass). 

Our results are given in section [5] and we quote from there a value for the 
nucleon sigma terms of 

a;"^)* = 31(3)(4)MeV, af > = 71(34)(59) MeV . (41) 

(The first error is the fit error while the second error indicates possible effects 
from higher order terms in the flavour expansion.) Note that expansions about 
the SU{3) flavour line require consistency between many QCD observables, here 
for example not only for the baryon octet under consideration here, but also for 
the pseudoscalar octet, and PCAC and the ratio of the light to strange quark 
mass. 
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Of course there are several more avenues to investigate. Numerically an in- 
crease in statistics for the masses along the flavour symmetric line would reduce 
the dominant error (both statistical and systematic) and so directly help in de- 
creasing the present errors. Our approach here has been to emphasise linearity 
at the expense (presently) of reaching exactly the physical point. This can be 
addressed by interpolating between a small set of constant m lines about the phys- 
ical point. Additionally the use of partial quenching will also help to get closer to 
the physical pion mass. With more data, a systematic investigation of quadratic 
quark mass terms in the flavour expansion should be considered, to reduce the 
systematic errors. Finally while the use of linear or quadratic terms along the line 
of constant m is unproblematic, so that it is unlikely that eq. ( HOj) will change by 
much, more subtle is the relation involving X{m) (i.e. the gradient when chang- 
ing m.) For the example of clover fermions we have g'^irn) = (1 + hgam)g'^ which 
clearly does not change if m = constant, but will slightly change when m does. 
However this is probably not a large effect (as hg seems small). For a discussion 
of some aspects of this issue see [291 EO] • 
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Appendix 

A Some relations between the a terms 

We discuss here some relations between the sigma terms within a multiplet (here 
taken to be the baryon octet) which are exact within the linear case discussed 
here, but which we might expect to always be approximately true. 

The singlet relation eq. fl22|) or eq. (j25|) is the same for every hadron. So in 
terms of sigma terms this becomes 

^\^) + raf^ ^ + raf^ . (42) 
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At the flavour symmetric point it follows from group theory that a singlet operator 
has the same value for every member of a multiplet, so eq. (142|) must hold. But 
this can change if we move away from the symmetric point. (We shall briefly 
discuss this at the end of this section.) 

We can find another collection of near identities by summing over a singlet 
combination of hadrons — this can be either a singlet of 6*3 or a singlet of SU (3). 
If we do this, the expectation values of uu, dd and ss will be exactly equal at the 
flavour symmetry point, and stay again nearly equal away from the symmetry 
point. By this argument we expect 

a^ + aP - 2r(af) + af)) 
-f' + <^'P + ^ - 2r (af ) + af ) + af )) . (43) 

(Again this relation, as with the other relations discussed here, is exactly true for 
the linear case.) 

Other relations come from the Gell-Mann-Okubo relation, [27J [28] in which 
the 27-plet mass combination is very small, 

2Mn - 3Ma - Me + 2Mh ^ , (44) 

for all values of mi,ms. In our approach, its derivatives are also near zero. We 
therefore expect 

2a(^) - 3a(^) - + 2aP ^ 

2af ) - 3a(^) - af ) + 2af ) ^ 0. (45) 

We obtain an even stronger version of these relations by taking the singlet com- 
bination, proportional to (uu + dd)'^ + (ss)^, 

2^;^) _ 3^;^) - ) + 2aP + r (2af ) - 3al^) - af ) + 2af )) ^ . (46) 

There is also a relation between the sigma terms and the hadron masses, 
[2] as the constants Ai and A2 which occur in the mass splittings also occur in 
the leading order expressions for the sigma terms. So there will be connections 
between masses and sigma terms. One particularly simple relation is 

Mh - of^ - o-f ) ^ Mh' - - ) . (47) 

(i.e. the baryon mass difference is closely accounted for by the sigma terms.) For 
the linear case this is again exact, with this equation being equal to Mo(mo) — 
moMQ(mo) for all the octet baryons (upon using eqs. ([7]), f l27|) ). From eq. f lTTj) we 
see that this is just the common hadron mass in the chiral limit along the flavour 
symmetric line, when m/ = = or m = 0. ct; and as can be thought of 
as that part of the hadron mass which is due to rrii and rris respectively. The 
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remnant, Mo(mo) — nioMQ^mo) , is the part of the hadron mass due to the quark 
and gluon kinetic energy, interaction energy, etc., [2], i.e. the part of the hadron 
mass which is not due to the couphng with the Higgs vacuum expectation value. 

We can use the higher order mass equations in [23j to estimate how well the 
relations in this section hold. Most of the relations have violations proportional 
to the first power of the SU{3) breaking parameter, Snii. The corrections to 
eqs. P2|) and fj43|) and the first relation in eq. fHSj) are 0{mi5mi). The as relation 
in eq. f H5|) has corrections 0{ms5mi). When we combine these two relations 
to form eq. f H6|) . the leading violation terms cancel, and we have a relation with 
corrections 0{miSmf). The corrections to the mass relation eq. ( H71) are 0{fnSmi) 
and 0{6mf). 

B Higher order effects 

In this Appendix, we discuss a little more quantitatively the systematic errors in- 
duced by the inclusion of the quadratic terms in the fit formulae. We concentrate 
particularly on the nucleon sigma terms, a^^^ and ai'^\ 

B.l Curvature in the 'fan' plot 

In Fig. |5]we compare the results of a quadratic fit and a linear fit, both for the 
baryon mass fan plot and on a[^^ and ai^\ The quadratic fit uses all the data, 
[23], on both lattice sizes (in cases where results for two lattice sizes are available, 
we used the larger lattice size only). Including curvature terms in eq. (^^, [23] , 
we have Mh = Mo + cnSmi + hn^Tn^ + . . .. Tracing through the analysis, we find 
the effect on eq. fl27|) is to replace 

ch ch + 2bH5mi . (48) 

By comparing Ch from the linear fit with Ch + 26//5m;* from the quadratic fit, we 
can estimate the maximum possible change. 

We use the data at kq = 0.12090, because this is the case where we have the 
most data, covering the largest range in quark mass splitting, 6mi. In this case 
we have data covering about 3/4 of the gap from the symmetric point to the 
physical point, so we have the most chance of seeing curvature effects if they are 
present. 

For the fan plot (left panel of Fig. [5]), the curvature terms are found to be 
small, and statistically compatible with zero curvature. In Fig. Owe compare the 
nucleon sigma terms from the slopes of the two fits by using eq. (^7^ together 
with eq. ( 148|) . Again we see that the curvature effect is very small in the case 
of particularly at small m^, and much larger for ai^\ Can we explain this 
difference? 
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Figure 9: (decreasing red lines from right to left) and Us (increasing blue lines 

from right to left) against M^/X^ using linear fits (dashed lines) and quadratic fits 
(solid lines) for kq = 0.12090. 



The slopes in the fan plot only effect the non-singlet matrix element, the ch 
term in eq. (127|) . The curvature changes the slope of the nucleon line by about 
10% at the physical point. The non-singlet term in a[^^ is responsible for about 
25% of the quantity, so a 10% change in slope translates to a 2.5% change in al^\ 
Putting in the actual slope change, the final number we arrive at is a systematic 
uncertainty of about 1 MeV in cr^ coming from curvature in the fan plot. 

The situation for ai^^ is different, the singlet and non-singlet terms appear 
with opposite signs, so ai^^ is given by the difference between two large quantities. 
Thus a 10% change in the non-singlet matrix element is leveraged into a 25% 
change in ai^\ Repeating this procedure for the other hadrons gives similar 
non-singlet uncertainties. 



B.2 Curvature along the symmetric line 

We also use a linear fit to describe the baryon masses along the symmetric line (the 
line with all three quark masses equal). What is the effect of using a quadratic 
fit to determine the slope along this line? 

In the right panel of Fig. |5] we compare a quadratic and linear fit to the 
symmetric baryon masses. As before, the quadratic term is compatible with zero 
curvature. Indeed the quadratic term is probably too large and is likely due to 
having a short lever arm and low statistics at the lightest point rather than to 
be a real effect. (Also we would expect that chiral perturbation theory would 
predict a downward curve.) 

Feeding these values into eq. ( |27|) gives an estimate of the possible effect of 
quadratic terms, due to curvature along the symmetric line, which we will include 
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in our final error estimate. This curvature effect is the same for every hadron, 
giving an uncertainty ~ 4MeV for ai and ~ 55MeV for a^. However because 
the shift is universal, this does not effect splittings, so the systematic error in 
al — cri is still given by the ~ 1 MeV value of the previous subsection. For 
y^^-*^, using the first equation in eq. (jlj) gives percentage changes in ?/(^)^ of 60% 
and 30% for and 

C Hadron Masses 

We collect here in Tables (E]) - ([5]) numerical values for the meson pseudoscalar 
octet and baryon octet, not given in [23j. (All the data sets used here are over 
~ 2000 configurations for the 24^ x 48 volumes and ~ 1500 — 2000 configurations 
for the 32^ x 64 volumes except for = 0.12099 which has ~ 500 configurations.) 
Errors are from a bootstrap analysis. 

~ 32^ X 64 

0.120920 0.1647(4) 0.4443(59) 



Table 2: Additional result for the pseudoscalar octet mesons and octet baryons along 
the flavour symmetric line: aMyr, aMjy, for (/3,Cs«,,a) = (5.50,2.65,0.1). 

{ki, Kg) gM^ gMx aM^^ 

24^ X 48 

(0.120870, 0.121020) 0.1804(8) 0.1621(10) 0.1407(12) 
(0.120980, 0.120800) 0.1545(9) 0.1775(8) 0.1976(7) 



Table 3: Additional results for the pseudoscalar octet mesons: aM^^, aMx and aM, 
for {P,Csw,a) = (5.50,2.65,0.1) where kq = 0.12092. 















24^ 


X 48 






(0.120870, 0.121020) 
(0.120980, 0.120800) 


0.4812(40) 
0.4668(61) 


0.4721(62) 
0.4773(62) 


0.4672(48) 
0.4838(47) 


0.4618(58) 
0.4909(41) 



Table 4: Additional results for the octet baryons: aMjv, oMa, aMs and aM= for 
{P,Cs^o,a) = (5.50,2.65,0.1) where kq = 0.12092. 
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[ku^^s) Mn/Xn Ma/Xjv M^/Xn Mh/X 



N 



X 48 


(0.120870, 0.121020) 
(0.120980, 0.120800) 


1.024(3) 1.004(9) 
0.9715(33) 0.9934(95) 


0.9939(17) 
1.007(2) 


0.9824(34) 
1.022(3) 


32^ X 64 


(0.121050, 0.120661) 


0.9167(40) 0.9872(46) 


1.017(2) 


1.066(3) 



Table 5: Additional ratio results for the octet baryons: M^/Xj^j, M\/Xj\f, My,/X^ 
and Ms/Xn for (/3,Cs^,a) = (5.50,2.65,0.1) where kq = 0.12092. 

References 

[1] T.-P. Cheng and L.-F. Li, Gauge Theory of Elementary Particles Oxford 
University Press (1988, reprinted); Schladming Winter School (March 1997), 
Computing Particle Properties (eds. C. B. Lang and H. Gausterer), Springer- 
Verlag, |arXiv : hep-ph/970 9293 1 . 

[2] X. Ji, Phys. Rev. Lett. 74 (1995) 1071, |arXiv : hep-ph/9410274] . 

[3] R. Koch, Z. Phys. C15 (1982) 161. 

[4] J. Gasser, H. Leutwyler and M. E. Sainio, Phys. Lett. B253 (1991) 252; 
ibid 260. 

[5] M. M. Pavan, L L Strakovsky, R. L. Workman and R. A. Arndt, PiN 
Newslett. 16 (2002) 110, | arXiv : hep-ph/0 111066 j . 

[6] B. Borasoy and U.-G. Meissner, Annals Phys. 254 (1997) 192, 
farXiv : hep-ph/9607432] . 

[7] J. Martin-Camalich, L. S. Geng and M. J. Vicente Vacas, Phys. Rev. D82 
(2010) 074504, [ arXiv : 1003 . 1929li rhep-lat1 ] . 

[8] J. Ellis, K. A. Olive and P. Sandick, New J. Phys. 11 (2009) 105015, 
jarXiv : 0905 . 0107i rhep-ph1 ] . 

[9] J. Giedt, A. W. Thomas and R. D. Young, Phys. Rev. Lett. 103 (2009) 
201802, [arXiv:090ai7T|[hep-ph]]. 

[10] S. Giisken, K. Schilling, R. Sommer, K. H. Miitter and A. Patel, Phys. Lett. 
B212 (1988) 216. 

[11] R. Altmeyer, M. Gockeler, R. Horsley, E. Laermann and G. Schier- 
holz, [MTc Collaboration], Nucl. Phys. Proc. Suppl. 34 (1994) 
376, arXiv : hep-lat/9311020[ Kyffhaeuser Workshop (Leipzig, Germany, 



September 1993), [iLrXrv:hep-lat/9311012| 



21 



[12] M. Fukugita, Y. Kuramashi, M. Okawa and A. Ukawa, Phys. Rev. D51 
(1995) 5319, [arXiv : hep-Tat/940 8002 1 . 

[13] S. J. Dong, J. F. Lagae and K. F. Liu, Phys. Rev. D54 (1996) 5496, 
[arXiv : hep-ph/9602259| . 

[14] S. Giisken, P. Ueberholz, J. Viehoff, N. Eicker, P. Lacock, T. Lippert, K. 
Schilling, A. Spitz and T. Struckmann, [SESAM Collaboration], Phys. Rev. 
D59 (1999) 054504, |arXiv : hep-lat/ 9809066 1 . 

[15] A. Walker-Loud, H.-W. Lin, K. Orginos, D. G. Richards, R. G. Edwards, 
M. Engelhardt, G. T. Flemming, Ph. Hagler, B. Musch, M. F. Lin, H. B. 
Meyer, J. W. Negele, A. V. Pochinsky, M. Procura, S. Syritsyn, C. J. Morn- 
ingstar, D. B. Renner and W. Schroers, Phys. Rev. D79 (2009) 054502, 
[a rXiv: 0806 .4549 [hep- 1 at] ] . 

[16] H. Ohki, H. Fukaya, S. Hashimoto, T. Kaneko, H. Matsufuru, J. Noaki, 
T. Onogi, E. Shintani and N. Yamada, [JLQCD Collaboration], Phys. Rev. 
D78 (2008) 054502, [arXiv : 806 . 4 744 fhep-lat! ] . 

[17] R. D. Young and A. W. Thomas, Phys. Rev. D81 (2010) 014503, 
[ arXiv: 0901 .3310 fhep-lat! ] . 

[18] K.-L Ishikawa, N. Ishizuka, T. Izubuchi, D. Kadoh, K. Kanaya, Y. Ku- 
ramashi, Y. Namekawa, M. Okawa, Y. Taniguchi, A. Ukawa, N. Ukita 
and T. Yoshie, [PACS-CS Collaboration], Phys. Rev. D80 (2009) 054502, 
|arX iv : 0905 . 0962 irhep-lat1 ] . 

[19] D. Toussaint and W. Freeman, [MILC Collaboration], Phys. Rev. Lett. 103 
(2009) 122002, [arXiv : 0905. 2432 [hep-lat]]. 

[20] K. Takeda, S. Aoki, S. Hashimoto, T. Kaneko, J. Noaki and T. 
Onogi, [JLQCD Collaboration], Phys. Rev. D83 (2011) 114506, 
[arXiv: 1011. 1964 [hep-lat] ] . 

[21] S. Diirr, Z. Fodor, T. Hemmert, C. Hoelbling, J. Prison, S. D. Katz, S. Krieg, 
T. Kurth, L. Lellouch, T. Lippert, A. Portelli, A. Ramos, A. Schafer and 
K. K. Szabo, larXiv : 1109 . 4265 fhep-lati . 

[22] R. D. Young and A. W. Thomas, (4th Int. Symposium on Symmetries in 
Subatomic Physics (SSP2009), Taipei, Taiwan, June 2-5 2009), Nucl. Phys. 
A844 (2010) 266C-271C. farXiv : 0911 .TTSTi rhep-latl . 

[23] W. Bietenholz, V. Bornyakov, M. Gockeler, R. Horsley, W. G. Lockhart, Y. 
Nakamura, H. Perlt, D. Pleiter, P. E. L. Rakow, G. Schierholz, A. Schiller, 
T. Streuer, H. Stiiben, F. Winter and J. M. Zanotti, [QCDSF-UKQCD Col- 
laboration], Phys. Rev. D84 (2011) 054509, [arX iv : 1 102 . 5300i fhep-lat1 ] . 



22 



[24] N. Cundy, M. Gockeler, R. Horsley, T. Kaltenbrunner, A. D. Kennedy, Y. 
Nakamura, H. Perlt, D. Pleiter, P. E. L. Rakow, A. Schafer, G. Schierholz, 
A. Schiller, H. Stiiben and J. M. Zanotti, [QCDSF-UKQCD Collaboration], 
Phys. Rev. D79 (2009) 094507, [a rXiv : 0901 . 33 02 fhep-latl ] . 

[25] M. Gockeler, R. Horsley, A. C. Irving, D. Pleiter, P. E. L. Rakow, G. Schier- 
holz and H. Stiiben, [QCDSF-UKQCD Collaboration], Phys. Lett. B639 
(2006) 307, [ arXiv : hep-ph/0409312j . 

[26] P. E. L. Rakow, Nucl. Phys. Proc. Suppl. 140 (2005) 34, 
larXiv-hep^lat/0411036] . 

[27] M. Gell-Mann, Phys. Rev. 125 (1962) 1067. 

[28] S. Okubo, Prog. Theor. Phys. 27 (1962) 949. 

[29] C. Michael, C. McNeile and D. Hepburn, Nucl. Phys. Proc. Suppl. 106 
(2002) 293, [ arXiv : hep-lat/0109028] . 

[30] R. Babich, R. C. Brower, M. A. Clark, G. T. Fleming, J. C. Osborn, C. 
Rebbi and D. Schaich, iarXiv: 1012.0562i 

[31] Y. Nakamura and H. Stiiben, 

PoS(Lattice 2010) 040, !!arMvil011^199Thep-lat] . 

[32] P. A. Boyle, Comp. Phys. Comm. 180 (2009) 2739. 

[33] R. G. Edwards and B. Joo, Nucl. Phys. Proc. Suppl. 140 (2005) 832 
[arXiv : hep-lat/0409003j . 



23 



